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Abstract 

We investigate the incompressible Navier-Stokes equations with variable density. 
The aim is to prove existence and uniqueness results in the case of discontinuous ini- 
tial density. In dimension n = 2, 3, assuming only that the initial density is bounded 
and bounded away from zero, and that the initial velocity is smooth enough, we get 
the local-in-time existence of unique solutions. Uniqueness holds in any dimension 
and for a wider class of velocity helds. Let us emphasize that all those results are 
true for piecewise constant densities with arbitrarily large jumps. Global results are 
established in dimension two if the density is close enough to a positive constant, and 
in n-dimension if, in addition, the initial velocity is small. The Lagrangian formula- 
tion for describing the flow plays a key role in the analysis that is proposed in the 
present paper. 
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Introduction 



Incompressible flows are often modeled by the homogeneous Navier-Stokes equations : 
that is the density of the fluid is assumed to be a constant. However in many applications as 
blood flows or models of rivers, although the fluid is practically incompressible, the density 
can not be considered as a constant quantity, as a consequence of the complex structure of 
the flow due to e.g. a mixture of fluids or pollution (see e.g. [2, 4, 19, 23, 29]). This makes 
us look at the density as a nonnegative unknown function which has constant values along 
the stream lines. The simplest model which can capture such a physical property is the 
so-called inhomogeneous Navier-Stokes system: 

Pt + v .^/p = in fix(0,T), 

pv t + pv ■ Vv - vAv + VQ = in ft x (0, T), 

divw = in fix(0,T), (0.1) 

on dn x (0,T), 



|t=o = v , p\ t =o = Po in Q. 



The unknown functions are: p - the density of the fluid, v - its velocity field and Q - its 
pressure. The constant positive viscosity coefficient is denoted by v. We consider the cases 



where Q is a bounded domain of IR n , or the whole space MJ 1 , and we focus mainly on the 
physically relevant space dimensions n = 2,3. 

The goal of the present paper is to revisit results concerning the well-posedness issue 
of (0.1). We concentrate our analysis on the regularity of density. In our recent work [11], 
we established the existence and uniqueness of solutions to (0.1) in a critical regularity 
framework which allowed the initial density to be discontinuous. However, a smallness 
condition over the jumps was required there. In the present work, we want to discard this 
smallness condition. At the same time, to simplify the presentation, we do not strive for 
optimal assumptions as concerns the velocity and assume the viscosity coefficient v to be 
constant. 

Let us recall (see in particular [13] and the textbooks [3, 20]) that, roughly, from the 
qualitative viewpoint the classical results for the homogeneous Navier-Stokes equations 
carry out to (0.1) : on the one hand global (possibly non unique) weak solutions with finite 
energy may be constructed and on the other hand, if the density is smooth enough, bounded 
and bounded away from zero, then global-in-time existence and uniqueness results are 
available in dimension two for arbitrarily large data, and if the velocity is small in dimension 
three. These latter results require relatively high regularity of the density, though. In 
particular it has to be at least continuous, and to have some fractional derivatives in 
suitable Lebesgue spaces (see e.g. [21] or [9]). It is worthwhile to emphasize that for 
smooth densities one may show the existence of unique solutions even for vacuum states 
[6] provided the initial data satisfy some compatibility condition. From the viewpoint 
of applications such results are not so satisfactory: we wish to consider fluids with e.g. 
piecewise constant densities, a pattern which is of interest to model a mixture of two 
fluids. 

The results of the paper are split into two groups: 

• The first group concerns uniqueness and local-in-time existence results in the case 
where the initial density is just an L^ positive function bounded away from zero. In 
particular, one may consider piecewise constant densities with arbitrary large jumps. 
As regards the existence issue, we have to restrict ourselves to the (physically relevant) 
dimensions n — 2, 3. 

• The second group concerns the global-in-time existence issue. Here we have to make a 
smallness assumption over the density which, in the case of piecewise constant initial 
density, implies that the jumps have to be small. Assuming enough smoothness over 
the velocity, this enables us to prove global existence for (possibly) large velocity if 
n — 2, and for small velocity if n > 3 (an assumption which is also required for the 
homogeneous Navier-Stokes equations, anyway). 

As explained above, in the present paper, we aim at doing minimal assumptions over the 
density but we do not strive for optimal regularity of the velocity function. As our method 
relies on estimates for the Stokes system with merely bounded coefficients, the (rather 
high) regularity of the velocity is somehow prescribed by the technique. An approach to 
the issue of sharp regularity has been done in [11] in the critical Besov spaces setting. 
However in [11] we were able to capture discontinuous density with small jumps only. 



The rest of the paper unfolds as follows. The main results are presented in the first 
section. Then, some preliminary estimates involving the evolutionary Stokes system are 
proved. Section 3 is devoted to the derivation of System (0.1) in Lagrangian coordinates. 
In Section 4, we concentrate on the proof of uniqueness results whereas existence results 
are proved in the last two sections. Technical estimates involving the divergence equation 
are presented in the Appendix. 

1 Main results 

Let us first recall the basic energy equality for System (0.1) which may be (formally) 
derived by testing (0.1)2 by v. 

Lemma 1 Let (p,v) be a sufficiently smooth solution to (0.1) over Q x [0, 7]. Then there 
holds 

(p\v\ 2 )(t,x)dx + 2u / \Vv(T,x)\ 2 dxdr = (p\v\ 2 )(0,x) dx for all t E [0,7]. (1.1) 
n Jo Jn Jn 

Subsequently if p is positive and bounded away from zero and v is in L 2 (£l) then we get 
a control over v in Loo(0, T; L 2 (fi)) and Vi> in L 2 (fi x (0, 7)). Under very rough regularity 
assumptions (much less than assumed here), the (formal) energy equality (1.1) provides 
us with an information about low norms of the velocity, which turns out to be crucial 
for the proof of global results (see in particular the monograph by [20] and the references 
therein, as regards the proof of global weak solutions with finite energy). Note that (1.1) 
gives some regularity information over the velocity even for very rough density. We shall 
see further in the paper a way to get even more regularity information over the velocity 
without assuming more on the density. 

Before listing the main results of the paper, let us introduce a few notation. Concerning 
the derivatives of functions / depending on both the time variable t and the space variable 
x, we denote by f t the time derivative and by Df the Jacobian matrix of / with respect 
to the space variable, namely (Df)ij = djf 1 . The notation V/ is reserved for T (Df). 

The Lebesgue spaces of measurable functions with integrable p-th power is denoted by 
L p (Q). More generally, if m E N then W™(Q) denotes the set of L p (Q) functions with 
derivatives of order less than or equal to m in L p (Q). Since the Navier-Stokes equations are 
of parabolic type, it is also natural to introduce parabolic Sobolev spaces W 2 ' p (Q x (0,T)) 
that is the closure of smooth functions for the norm 

H U llw, a ,' p 1 (nx(0,T)) = H U > d M\L p (0,T;L q (Cl)) + \\u\\ L p (0,T;W^(Q))- (1-2) 

Granted with parabolic spaces, one may now define Besov spaces over Q as the following 
trace space: 

B 2 ~ 2/p (n) = [/ : n ->■ R measurable s.t. / = f\ t=0 for some feW^(il x (0, 1))|- (1.3) 

The norm can be defined from the above definition as a suitable infimum (for more details 
concerning the Besov spaces we refer to [5, 33]). 



Our first result states the uniqueness of solutions with merely bounded density, provided 
the initial velocity is smooth enough. 

Theorem 1 Let n > 2. Assume that Q is M. n or a C 2 bounded domain of MJ 1 . Let 
(p 1 ,v 1 ,(5 1 ) and (p 2 ,v 2 ,Q 2 ) be two solutions to (0.1) with the same initial data, and density 
bounded and bounded away from 0. Suppose moreover that for k = 1,2, 

• Case n = 2: there exists q > 2 such that v k G W^'^Vt x (0,T)) and VQ k G 
L 2 (0,T;L q (n)), 

• Case n > 3: v k G W*i(Q x (0,T)) ; VQ k G L 2 (0,T; L n (fi)) and, in addition, Vv k G 
L 2 (0,T;L oo (n)). 

Then v 1 = v 2 , VQ 1 = VQ 2 and p 1 = p 2 . 

Remark 1 As regards the inhomogeneous incompressible Navier-Stokes equations, to our 
knowledge, the "best" uniqueness result with no smallness condition over the density is due 
to P. Germain in [16]. It does not apply to solutions with piecewise constant densities, 
though. 

The second result complements Theorem 1. It delivers existence of local-in-time regular 
and unique solutions in dimensions 2 and 3. Again, the initial density just has to be bounded 
and bounded away from vacuum. 

Theorem 2 Let n = 2, 3 and Q be a C 2 bounded domain or be W 1 . Let p satisfy 

m< p < M (1.4) 

for some positive constants, and Vq G W^ity be such that divt> — and t>o|an = 0. Let 
n* = 2( ! ^). There exists a unique solution (p,v) to System (0.1) on a time interval [0, T] 
for some T > such that p(t, •) satisfies (1.4) for all t G [0, T] and 

veW$ n .(nx(0,T)), v t G Loo^T; L 2 (Q)) and Vv t G L 2 (fi x (0,T)). 

Remark 2 The critical Sobolev embedding ensures that W$(Q) is continuously embedded 
in the Besov space B n 7„* (fi). Keeping in mind the definition of this space given in (1.3), 
the appearance of the parabolic Sobolev space W n l n * (Q x (0, T)) in the above statement 
does not come up as a surprise. The W 7 ! (fi) assumption for vq is needed to ensure that 
(dtv + v ■ Vi>)|t=o is in L 2 (fi). At the same time, owing to the low regularity of the density, 
we do not know how to propagate the W 2 (Q) regularity for the velocity. 

Proofs of Theorems 1 and 2 are based on the analysis of (0.1) in the Lagrangian coordi- 
nates defined by the stream lines. Since the density is merely bounded there is an obstacle 
to apply any bootstrap method in order to improve the regularity of the velocity. The 
main difficulty is located in the term with the time derivative. To obtain a better infor- 
mation about v t , we adopt techniques from the compressible Navier-Stokes system [25, 28] 
(concerning uniqueness criteria for the compressible Navier-Stokes system in Lagrangian 
formulation, see also the recent work by D. Hoff). Roughly speaking, we differentiate the 



(Lagrangian) velocity equation once with respect to time, then apply an energy method. 
This approach via the Lagrangian coordinates requires only L^ bounds (by above and 
by below) for the density, provided the velocity has high regularity. That the density is 
time-independent in the Lagrangian setting, hence is just a given function, is of course fun- 
damental. In comparison, in [9] where the Eulerian framework is used, the initial density 
has to be in the Besov space B^ x (W n ) (which, roughly, means that it has n/2 deriva- 
tives in L 2 (M. n )) but the initial velocity therein has only critical regularity, namely it is in 
i?2 i ~ (K n ) (to be compared with B 2 1 (M n ) and n = 2, 3 here). 

To highlight consequences of Theorems 1 and 2, let us consider the case where the 
initial divergence- free velocity field is in W%(VL) (and vanishes at the boundary), and the 
initial density p u is 1 

Po =m + axA , (1.5) 

where m, a are positive constants and A$ is a set with a C 1 boundary. The velocity field 
v given by Theorem 2 is Lipschitz with respect to the space variable hence generates a 
unique C 1 flow X defined by 

X(t,y)=y+ [ v(r,X(r,y))dr. 
Jo 

Therefore, the density at time t is given by 

p(t,-) = m + a X A(t), with A(t):=X{t,A ). (1.6) 

As the flow X is at least C 1 , the initial regularity of the boundary of A(t) is preserved and 
any geometrical catastrophe (e.g. breaking down or self-intersections of the boundary) will 
not appear : if Aq is diffeomorphic to a ball, then A(t) is diffeomorphic to a ball, too. The 
above case shows that the system (0.1) can model an interaction of two fluids separated by 
a free interface. Although tracking the regularity of the boundary dA(t) is not the main 
topic of this paper, we see that Theorem 2 ensures that the C 1 or C 1,a regularity (with a 
small enough) of dA(t) is preserved during the evolution. In other words, we have partially 
solved in an indirect way a complex free boundary problem which has been left as an open 
question by P.-L. Lions in his book [20]. Let us emphasize that the standard approach for 
solving problems of such type requires very technical considerations (see e.g. [1, 30, 32]). 
Furthermore, with our approach, there are no requirements concerning the regularity of 
the boundary of the set Aq: our results hold for any measurable set Aq. 

The above results concern local-in-time analysis. In order to obtain global-in-time so- 
lutions, we have to assume that the jumps of the initial density are small enough. The 
following theorem states that under this sole assumption over the density, and for suffi- 
ciently smooth (possibly large) initial velocity fields, global existence holds true. 

Theorem 3 Let Q be a C 2 bounded two-dimensional set, or be M?. There exists a constant 
c depending only on Q and such that if po G L^IQ) satisfies 

suppo(^) — inf Po(x) 
e —. n < c (1.7) 



^cre xa stands for the characteristic function of the set Aq. 



then for all vq e B\ 2 {&) fl L 2 {0) with divt>o = and t>o|an = 0, there exists a unique 
global-in-time solution to System (0.1) such that (1.1) is satisfied and that, for allT > 0, 

veW^ 2 (Qx(0,T)), VgeL 2 (0,T;L 4 (fi)) and p e L^Q x (0,T)). (1.8) 

In dimension n > 3, getting global-in-time strong solutions requires also the initial 
velocity to be small, (an assumption which is needed in the homogeneous case, anyway). 
Here is our statement: 

Theorem 4 Letfi be a bounded n- dimensional C 2 domain. Let p e L 0O (Q) be positive and 

bounded away from 0, and v G -£>q, P 9 (fi) with 1 < p < oo, n < q < oo and 2 — 2/p ^ 1/q. 
There exist two constants c and d depending only on Q, p and q and such that if 

suppo(^) — hif po(x) 

X€z£l Xdzll t II II / / \ 

<c and \\vq\\ 2 _2 <cu, (1.9) 

ml p (x) B q , p p (n) 

then there exists a unique global-in-time solution to the inhomogeneous Navier-Stokes sys- 
tem (0.1) such that 

v e W^(n x R + ), VQ e L P (R + ; L q (Q)) and p e L^Q X R+). 

Furthermore, there exist two positive constants a and C depending only on Q, p, q and of 
the lower and upper bounds for po so that for all t > 0, 

W v \\w*£(cix(t,t+i)) + l|VP||L p (t,t+i;L,(n)) < Ce~ a *||f;o|| 2 _| . 

Theorems 3 and 4 follow from classical maximal regularity techniques. The smallness 
conditions (1.7) and (1.9) allow to treat the oscillations of the density as a perturbation 
that may be put in the right-hand side of the estimates. 

At the end we would like to underline that most of our results hold for bounded domains 
and M™. The case of the whole space is easier: there is no boundary condition and solving 
the divergence equation is simpler, too. One exception is Theorem 4 where the boundedness 
of the domain is essential here as it provides exponential decay of the energy norm (the 
whole space case is tractable under stronger conditions over the density, see our recent 
work in [11]). 

2 Some linear estimates 

A fundamental role will be played by the Stokes system, that is the following lineariza- 
tion of the velocity equation in (0.1): 

mu t -vAu + VQ = f in fix(0,T), 

divM = divi? in Qx(0,T), 

V ; (2.1) 

u = in fflx (0,T), 

u\ t =o = Mo in fi, 

where m and v are positive constants. 



We shall make an extensive use of the following solvability result for the Stokes system 
in the L p (0,T; L q (Q)) framework: 

2-2 

Theorem 5 Let Q be a C 2 bounded subset ofW 1 . Let 1 < p,q < oo, u <E B QtP p (fi), 
f G L p (0,T;L q (tt)), R G W}(0,T;L q (Q)) so that div R G L p (0,T;W}(Q)). Suppose that 

divu = div-R| i= o and ft- R\onx(o,T) = 0. 

If 2 — 2 /p > 1/q, assume in addition thatuo = at the boundary, otherwise we assume only 
u -n = at dQ 2 . Then there exists a unique solution to (2.1) such thatu G Wg£(Qx(0, T)), 
VP G L p (0,T; L q (Q)), and the following estimate is valid: 



\\mut, z/V u, VP||r (o,T-,L q (n)) +mpi/p' sup ||u(^n 2 _, 

0<t<T B qiP p (Q) 

l i 

<C(\\f,mR t \\ Lpi0tT ; Lqin)) + \\vdwR\\ Lpi0:T . w i m +mpv7\\u \\ 2 _a ), (2.2) 

Bq,p (f2) 

where C is independent of m, v and T . 

Proof: In the case R = 0, this result is classical (see e.g. [17, 24] and the appendix of [7]). 
The general case follows from this particular case once constructed a suitable vector-field 
w : il x (0, T) ->■ W 1 fulfilling 

divw = divi? in fi, ty = at dQ. (2.3) 

Taking for granted the existence of such a vector-field, the strategy is simple : we set 
v = u — w and we gather that v has to satisfy 

mv t — uAv + VQ — f — mw t + uAw in fix (0, T), 

divw = in fix(0,T), 

w = in dttx (0,T), 

v|t=o = Mo-w in fi. 

Therefore, in order to reduce our study to the case R = 0, the vector-field u> is required 
to satisfy w t , D 2 w G L p (0,T; L q (£l)) (note that this will imply that w G P g ,p (fi), see 
(1.3)). The fact that such a solution to (2.3) does exist is granted by the following lemma, 
the proof of which is postponed in Appendix (see Proposition 3): ■ 

Lemma 2 Let R(t, ■) be a family of vector-fields defined over the C 2 bounded domain Q, 
parameterized by t G (0,T). Assume that, for some 1 < q < oo and 1 < p < oo we have 
div Re L p (0,T;W*(n)), R,R t G L g (0,T;L p (fi)) and R ■ n = at the boundary. 
Then there exists a vector-field w in L p (0,T; W 2 (Q)) vanishing on dQ, fulfilling 

divw = div R and divu; t = divi? t in Q 



2 For simplicity we exclude the case 2 — 2/p = 1/q. 



and the following estimates: 

||w||Lp(o,T;Wj(ft)) < C\\divR\\ Lp{0yT . w i m , (2.4) 

\\Wt\\L p (0,T;L q (n)) < C\\R t \\L p (0,T;L q (n)) (2.5) 

for some constant C depending only on q and Q. 

Remark 3 The whole space case is easier to deal with for we do not have to take care of 
boundary conditions (apart from suitable decay at infinity given by the functional setting). 
Indeed, in order to solve (2.3), one may set 



w 



-V(-A) _1 divP. 



As the corresponding Fourier multiplier is homogenous of degree 0, we readily get (2.4) 
and (2.5). Therefore, arguing as above and using the standard maximal regularity result 
for the Stokes system in IR n , we conclude to Theorem 5 in the case Vt = W 1 if the Besov 



2-- • 2- - 

space B q<p p (Q) is replaced with the homogeneous Besov space B qtP v (W 1 ) and W£(Q), by 
its homogeneous version Wg(W l ). 

Theorem 5 can be viewed as a classical result. In order to prove Theorem 2 we need to 
adapt it to the case of variable coefficients. Below, we focus on the L 2 case where only the 
boundedness of coefficients is needed. 

Lemma 3 Let Q be a bounded domain ofM. n , orM. n . Let r\ G L 0O (fi) be a time independent 
positive function, bounded away from zero, and R satisfy the above boundary conditions. 
Then the solution (u, VP) with u G W^{Q x (0,T)) and VP G L 2 (fi x (0,T)) to the 
system 

r]u t - vAu + VP = / in fix(0,T), 

divu = divP in Q x (0, T), , , 

u = on <9fix(0,T), *• ' 

u\t=o = Mo on Vt, 

fulfills 

^/v SUp \\Vu(t)\\ L2{ n) + ||M t ,Z/V 2 M,VP|| L2 (Qx(0,T)) < C(||/,P<||l 2 (Qx(0,T)) 
0<<<T 

+ v\\<^ R\\l 2 {o,t-w^)) + v^||Vuo||L a (n)), (2.7) 

where C depends on infi] and sup 77, but is independent ofT and v. 

Proof: First we remove the right-hand side of (2.6) 2 by means of Lemma 2 (or the 
remark that follows if Q = M. n ): we introduce a vector-field w fulfilling (2.3) such that 
w G W 2 ' 2 (^ x (0, T)) with the following bound 

\\dtw, vD 2 w\\ L2{ o,T;L 2 (n)) < C(\\dwR\\ L2i0:T . w i {n)) + \\R t \\L 2 (nx(o,T)))- (2.8) 

Hence we may reduce the proof to the case R = 0. Now, we observe that testing by u t 
gives: 

, W ^ + ^/|v u |^=//.M*. 

q 2 at J n J n 



Therefore, integrating in time yields 

"\\Vu(t)\\l m + / WVvutWl^dr < v\\Vu \\l 2{n) + [ \\f/Vv\\l 2(n) dr. (2.9) 

Jo Jo 

Since rj is a positive time independent function which is pointwise bounded from below 
and above, we obtain 

IMk(nx(o,T)) + sup y/v\\Du(t)\\ L2 (n) < C(||/lk(£}x(o,T)) + y/u\\Du \\ L ^n)). (2.10) 

0<t<T 

In order to estimate D 2 u and DP, we rewrite (2.6) as 

-uAu + VP = f -r]u t in fix(0,T), 

divw = in ftx(0,T), (2.11) 

m = at <9fix(0,T). 

If n is a C 2 bounded domain then the solvability of (2.11) in the L 2 framework is 
clear (see e.g. [15], Th. 6.1, page 231), thus taking into account bounds (2.8) and (2.10) 
we get (2.7). Lemma 3 is proved. In the IR n case, one may just notice that VP = 
- V(- A) _1 div (/-??«*). As f-u t is in L 2 {Vt x (0,T)), we still get the result, first for VP, 
and next for V 2 w. ■ 



3 The Lagrangian coordinates 

A fundamental point of our analysis is the use of Lagrangian coordinates. In order to 
define them we solve the following ordinary differential equation (treating y as a parameter): 

^M = w (t,X(t,y)), X(t,y)\ t=0 = y. (3.1) 

This leads to the following relation between the Eulerian coordinates x and the Lagrangian 
coordinates y: 

X(t,y) = y+ [ v(r,X(r,y))dr. (3.2) 

Jo 

Let us list a few basic properties for the Lagrangian change of variables: 

Proposition 1 Suppose that v G 1^(0, T; W^Q)) with v ■ n|g^ = 0. Then the solution 
to System (3.1) exists on the time interval [0, T], X{t,Vt) = Q for all t G [0,T), and 
D y X G Loo(0,T; L 0O (fi)) with in addition 

\\D y X(t)\\ Loc{ n)<exp( \\D x v\\ Loc(n) drj . (3.3) 



Furthermore 



X(t, y)=y+ [ u(t f , y) dt! with u(t, y) := v(t, X(t, y)) (3.4) 

Jo 



so that DX satisfies 



(3.5) 



D v X(t, y) = Id + / D y u(t',y)dt'. 
Jo 

Let Y(t, ■) be the inverse diffeomorphism of X(t, •). Then 

D x Y(t,x) = (D y X(t,y))- 1 with x = X(t,y) (3.6) 

and, if / \D y u(t,y)\ dt! < 1/2 then 
Jo 

\D x Y(t, x) - Id| < 2 / |Z? y w(i', y) | dt'. 
Fma% 7 i/w G Li(0,T; W^(Q)) then D y X G Loo(0,T; W^(Q)) and 

|^X(t,y)|< e /ol^(''^'^))l d '' f\Dlv{t\X{t',y))\eti\ D *< t "> x «'^\ dt "dt'. (3. 



(3.7) 



andi/v e Li(0,T; W£(ft)) with s > J + 1, tfien Dj,X - Id G £«,((), T; W^ _1 (fi)). 

Proof: The existence of X for (£, y) G (0,T) x f2 follows from the standard ODE theory, 
a consequence of Picard's theorem. Inequalities (3.3) and (3.8) follow from (3.2) by dif- 
ferentiation and Gronwall lemma. The higher regularity stems from the fact that, under 
our assumptions, W^ 1 is an algebra (the reader may refer to the appendix of [11] for the 
proof of similar results in a slightly different context). 

Equation (3.5) follows from (3.4), by differentiation. Then (3.7) comes from (3.6) 
provided D y X — Id is small enough: indeed, we have 

+oo , „ t -, fc 

D X Y = (Id +(D y X -Id))- 1 = J2(- l ) k [ D y u(t',y)dt'\ . 

This yields (3.5). ■ 

Let us now derive the Navier-Stokes equations (0.1) in the Lagrangian coordinates: we 
set 

r)(t,y):=p(t,X(t,y)), u(t,y) := v(t,X(t,y)) and P(t,y) := Q(t,X(t,y)). (3.9) 

We claim that System (0.1) recasts in 

r] t = in fix(0,T), 

r]u t - uA u u + V U P = in fix(0,T), 

div u w = in fix(0,T), (3.10) 

u = on dtt x (0,T), 

u\ t =t = v\ t=to , i]\ t=0 = p\ t=to in Q, 

where operators A u , V u , div u correspond to the original operators A, V, div , respectively, 
after performing the change to the Lagrangian coordinates. Index u underlines the depen- 
dence with respect to u. Let us also notice that, as v and u vanish at the boundary, we do 
have X(t,Q) = D, for all t. 

10 



So let us now give a formal derivation of (3.10). First, given the definition of X, it is 
obvious from the chain rule that 

d t r)(t, y) — (d t p + v ■ Vp)(t, x) and d t u(t, y) — (d t v + v • Vv)(t,x) with x :— X(t,y). 

The chain rule also yields 

D y P(t,y) = D x Q(X(t,y))-D y X(t,y) with (D y X) tj := d y .X\ (3.11) 

Hence we have 

D x Q(t, x) = D y P{t, y) ■ A(t, y) with A(t, y) := (D y X(t, y))- 1 = D x Y{t, x). (3.12) 

Next, we notice that if the transform X is volume preserving then for any smooth 
enough vector-field H we have 

div x H(x) = dW y (AH)(y) with x = X(y) and H(y) = H(x). (3.13) 

This stems from the following series of computations which uses the fact that det A = 1 
and the change of variable x = X(y): for any smooth q with compact support, we have 

/ q(x) div,,, H(x) dx = — D x q(x) ■ H(x) dx, 

D y q{y) ■ A{y) ■ H{y) dy, 

q(y)div y (AH)(y)dy. 



Combining (3.12) and (3.13), we thus deduce that, in Lagrangian coordinates operators 
V, div and A become 

V U := T A-V y , div u := div (A-) and A u :=div u V u . (3.14) 

In consequence, we have the following relations that will be of constant use: 

(V - V U )P = (Id - T A) VP, (3.15) 

(A-A u )w = div((Id-A T A)VM). (3.16) 

Let us finally emphasize that, owing to the chain rule, we have 

dw y (A-) = div M = A: D y . (3.17) 

This algebraic relation will be of fundamental importance in our analysis. 

The following statement ensures the full equivalence between (0.1) and (3.10) under 
the assumptions of our results stated in Section 1. 

Proposition 2 Let 1 < p, q < oo. Let p G Loo(^) and (u,P) be a solution to (3.10) such 
that u E W^(Q x (0,T)), VP G L p (0,T;L q (Q)) and 

-T 

II Vu|| ioo( n) dt < 1/2. (3.18) 



Then 

v(t, x) = u(t, y), Q(t, x) = P(t, y) and p(t, x) = p (y) 

with x = X(t,y) given by (3.2) defines a W^' 1 -solution to (0.1). 
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Conversely, if p E L^Q x (0,T)) and (v,Q) with v G W 2 £(Q x (0,T)), Vv e 
Li(0,T; L 0O (Q)), and VQ G L p (0,T; L q (Q)) is a solution to (0.1) t/ien 

u(t, y) = v(t, X(t, y)), P(t, y) = Q(t, X(t, y)) and rj = p\ t=0 

defines a W 2, p -solution to (3.10). 

Proof: The proof goes along the lines of the corresponding one in the appendix of [1 1] . 
Having Du small enough in Li(0, T; L 0O (fi)) is of course fundamental. ■ 

4 Proof of Theorem 1 — uniqueness 

In this part we prove the uniqueness of solutions to System (3.10) under the assumptions 
of Theorem 1. Here Q is a C 2 bounded domain, or the whole space. The proof is a 
straightforward application of Lemma 3 to the equations in the Lagrangian form. The 
important fact is that we have 

Wel^T;!^)) 2 = 1,2. 

Hence, taking T small enough, one may assume with no loss of generality that 



|Vu*|U oo( n ) dt< ^> 



(4.1) 



so that Propositions 1 and 2 apply. In particular the regularity properties of those solutions 
in Lagrangian coordinates are the same as those of Theorem 1. Hence it suffices to consider 
two solutions (u l , P 1 ) and (u 2 , P 2 ) to System (3.10) with the same initial data and satisfying 
the conditions of Theorem 1. 

Then, denoting A 1 := A{u l ) (see (3.12)), du := u l — u 2 and so on, we get 



r]5a t - z/A&z + V5P 



div&z = div [(Id 
6u = 
Su\ t= o = 



-z/[(A - A u iX - (A - A u 2)u 2 } 
+ [(V-V u i)P 1 -(V-V u2 )P 2 ] 

A^u 1 - (Id - A 2 )u 2 } 



in fix(0,T), 

in nx(0,T), 

at dQx(0,T), 

in Q. 



(4.2) 



Let us underline that the boundary condition on R from Lemma 3 is fulfilled, since by 
definition u l and u 2 are zero at the boundary. Therefore, keeping (3.17) in mind, we 
obtain for some constant C depending only on v, infrj, sup 77 and Q the inequality 



INUooto.TW^)) + ll&t, V 2 <5u, V5P|| L2( n X (o,T)) < C(h + I 2 + h + h) 



(4.3) 



with 



h 
h 
h 
h 



V u 2)P 2 || L2 ( nx ( 0iT )), 



IKv-v^p 1 -^ 

||(Id - A 1 ) : Du' - (Id - A 2 ) : Du 2 \\ W]W , m , 
|| [(A- A u iX - (A- A u 2)M 2 || L2 (n x( o,T)), 
||9 t [(Id - A l )u l - (Id - A 2 ) M 2 ]|| L2(nx(0 ,T)). 
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In the following computations, we shall use repeatedly the fact that SA := A 2 — A 1 
satisfies 



8A(t) 



[ Ddudr) ■ (j2 J2 CiCj" 1_i> ) with d(t) ■= [ 



i(t) := I Du'dr. (4.4) 

o 



We concentrate on the case n > 3. We shall indicate how our arguments have to be 
modified if n — 2, at the end of the section. 
In order to bound I±, we write 

/ 1 (i)<||(A 1 -A 2 )VP 1 |U 2( nx(o )i) ) + ||(Id-A 2 )V(P 1 -P 2 )|U 2( nx(o,t)). (4.5) 

It is clear that 

||(Id- A^ViP 1 -P 2 )\\ L2{ a Am < C* 1/2 ||VdP|| La( nx ( o,t))||i?ti 2 ||L a (o,t;L«(n)). (4.6) 

Let us notice that, according to (3. 7), (4. 4) and to the critical Sobolev embedding of W 2 (Q) 
in L 2 *(Q) (that is 1/2* + 1/ra = 1/2), we have 



¥A\\ 



Loo(0,t;L 2 *(n)) 



< CI 



\V8u\dt'\ 



Loo(0,t;L 2 *(n)) 



< Ct 1/2 \\D 2 ou\\ L2{Qx{ o,t)) 
with C depending only on the norm of the two solutions on [0, T\. Therefore, 

IKA 1 - A 2 )VP 1 \\ L2 (nx(o,t)) < C , ||&4|| Loo (o,t;L 2 »(n))||VP 1 ||L 2 (o,t;L n (n)) 

< Ct 1 / 2 ||D 2 Ai||L 2 (nx(o,t))||VP 1 ||L a (o ) t;L„(n))- 

Let us now bound I 2 . Note that it suffices to bound the norm in L 2 (Q x (0, T)) of the 
gradient of the corresponding term. If Q is bounded this is a consequence of the Poincare- 
Wirtinger inequality as div Su has average over Q, and if Q = W 1 this stems from the fact 
that only the norm in W^R") is involved (see Remark 3). Now, we notice that 

(Id - A 1 ) : Du 1 - (Id - A 2 ) : Du 2 = -5A : Du 1 + (A 2 - Id) : D6u. 
First, using the embedding of W 2 (Q) in L 2 *(Tl), and keeping in mind (4.1) and that 

Dtf e I 2 (0, T; ^(Q)) and DV £ L 2 (0, T; L n (Q)) fori = 1,2, (4.7) 

we get for all t e [0,7], 



|P>(&4: J D U 1 )|U 2( fi x(0 ,t))< \Du l \f \D 2 5u\dr' 

Jo 



L 2 (Clx(0,t)) 



+ 



P>V 



\D5u\ dr' 



L 2 (f2x(0,t)) 



< t 1/2 [\\D 2 5a\\ L , 2 (n x ^ t ))\\Du l \\ L2 (o,t;Loo(n)) 

+ ||-D(5u||L 2 (o,t;L 2 »(n))||-D 2 '" 1 ||L2(0,i;L„(n))J 
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Second, we have 
\\D((A 2 - Id) :£>«5u)|| La (nx(o ) t)) £ \\DA 2 ^D6u\\ L2(n>< ^ T)) + \\(A 2 - Id)®D 2 5u)\\ L2(n><m) 

< +1/2 ( \\ n2„,2|i ii n&.u 

~ l \\\ U u \\L2{0,t;L n (n))\\-L>dU\\ L2 (o tt -,L 2 *(n)) 

+ \\ Du2 \\L 1 (o,t;Lo C (n))\\ D2 fa\\L 2 (nx(o,t))y 
So finally for all £ G [0,T], 

/ 2 (t) < C* 1/2 ||V 2 «5u|U 2 (nx(o lt )) (4.8) 

with C depending only the norm of the solutions over [0, T\. The term J 3 may be handled 
along the same lines. Indeed we have 



/ 3 (t) = lldiv ((Id - A lT A l )Vu l - (Id - A 2T A 2 )Vu' 



L 2 (nx(o,t)) 



Finally, we examine J 4 . Using again (4.4), we get (with the convention that Du 1,2 denotes 
the components of Du 1 and Du 2 ): 



WdtlSAu^WL^nx^t))^ \\D5uu 1 \\ L2 ( nx ( 0tt ))+ / \D6u\dr'\Du ' \\u 

'o 



+ 



\D5a\dr'\u\\ 



L 2 (Qx(0,t)) 

L 2 (nx(o,t)) 



< H- D(5u IUoo(0,t;L2(^))ll Ml |U2(0,t;L oo (n)) 

+t 1 / 2 ||Z^^x|| jL2(0 ,t;^ 2 *(^)) (||^^ |U 2 (0,t ;J C„(^)) + || l^ 1 ! IZ^-U 1 ' 2 ! || jL2(0 ,t;i T .(^))) 

< t 1/2 ||D 2 &||L 2 (fix(0,i)) + £ WP^IUo O (0,t;L 2 (Q)), 

with lim t _ ) .o £ (^) — because 

SjM 1 , w^WandM^Dw 1 are in L 2 (0,T; L„(Q)). (4.9) 

At the same time, we have, for all t G [0,T], 

||<9 t ((Id - A 2 )6u)\\ L2{nx{ o :t)) < \\Du 2 Su\\ L2 (ax(o,t)) + ||(Id - A 2 )d t 5u\\ L2{Q ,x(o,t)) 

< INb 



|Loo(0,t;L n *(fi))||-Dw ||z, 2 (0,t;L„(f})) 



LDm 2 



iLi^LooCfi)) 



|W| 



L 2 (nx(o,t))- 



So one may conclude that 

hit) < t l/2 \\D 2 5u\\ L2{nxm) +e(t)(||-D&t||L 00 (o,t;L a (n)) + ||<9 t <5u||z, 2 (f2x(o,t)))- 

So finally in the case n > 3, putting together all the previous inequalities yields for all 
£G(0,T), ' 

INIwo^ffi)) + ll^t. V 2 *x, V5P|| L2(f7x(0 ,i)) 



< e(t)^||&i|| Loo(0i t ; wi(n)) + ll&*t, V 5u,V5P\\ L2{ nx(o,t)) 

for some positive function e going to at 0. Uniqueness follows on a sufficiently small time 
interval, then on the whole interval [0,T] thanks to a standard connectivity (or bootstrap) 
argument. 
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Let us now explain how the arguments have to be modified in the two-dimensional case. 
One cannot follow exactly the above approach owing to the failure of the embedding of 
W 2 (Q) in Loo(r2). So we have to assume slightly higher regularity, namely VP l ,VP 2 G 
L 2 (0,T; L q (Q)) with q > 2, and so on. For instance, setting m G (2, oo) such that l/m + 
1/q = 1/2, we may write 



IK^-A^VPi^oxfct)) < Ct^\\V5u\\L 2{ o,f,L mm \\VP* 



\L 2 {0,t;L q (Q,)) 



< C^2|| jD 2 (ii || i . 2(n><(0ii)) || V p2|| L2(o ^ (n)) _ 

The other terms of (4.3) may be handled similarly. The details are left to the reader. 
Theorem 1 is thus proved. 

Remark 4 Here we would like to explain the reason why we use the W 22 regularity for the 
velocity to establish uniqueness. Concentrate our attention on n = 3. A direct L 2 -energy 
method (i.e. testing (4.2) by 6u) requires our bounding (V — V^P 1 — (V — V U 2)P 2 in 
Li(0,T; L 2 {VL)), hence the following computation: 



MVP 1 Su dx dt 



< C||^4||L oo (0,T;L 2 (n))||VP 1 ||L 2 (0,T;L 3 (^))ll* i IU2(0,T;L 6 (Q)) 

< CT 1 / 2 ||VP 1 ||L 2 (0,T;L 3 (n))||V*i||L 2 (nx(0,T))- 



So we need VP 1 G L 2 (0,T; L 3 (f2)) which is naturally related to u l G W 32 . In addition 
integrating by parts in the left-hand side of the above inequality, we need to keep track of 
V 2 Su as well as ofV5P in L 2 (Q x (0,T)). Those two terms are out of control if resorting 
only to the basic energy inequality. 

5 Proof of Theorem 2 — existence 

The uniqueness property of the system is important, but to have the full picture of the 
well-posedness issue, we now have to show that there exist solutions with merely bounded 
density for which Theorem 1 applies. With the method that is proposed below, much more 
regularity is needed for the initial velocity. However the assumption over the initial density 
stays that same: it just has to be bounded and bounded away from zero. 

5.1 A priori estimates 

We first concentrate on the proof of a priori estimates for a smooth solution (u, P) to 
(3.10). To simplify the presentation, we consider the case where Q is a C 2 bounded domain 
of M. n . The whole space case may be achieved by similar arguments : this is just a matter 
of using homogeneous norms || ■ ||^i(ir™) an d || • \\^2-2/ n * ( „ n] and resorting to Remark 3. 

In order to prove a priori estimates for (u,P), let us assume in addition that T has 
been chosen so that (say) 



/ \\^u\\ Loo( n)dt<l/2. (5.r 
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This enables us to go from (0.1) to (3.10) (and conversely). For any (possibly large) initial 
velocity v e P n 7 n * (£1), and po £ ^oo(fi) bounded away from zero, we want to find a 
bound for a solution (u, P) given by Theorem 2. In other words, we want to control the 
following quantity: 

E(„,P)(T) := IKIUoo^TjLaCfJ)) + ||Vw t ||L 2 (nx(0,T)) + IK, V U, VP||x, n »(fi >< (o,T)), (5.2) 

with n* = 2( n ±^), if T is small enough. 

Let us first notice that, by standard Sobolev embedding 

/ HVulU^n) dt < CT 1 -^ Z {U , P) (T). (5.3) 

Jo 

which guarantees (5.1) for small times. 

In order to use Lemma 3 we restate System (3.10) as follows (of course t] = p and 
M = v ): 

rju t - vAu + VP = -i/(A - A u )u + (V - V U )P in fix(0,T), 

divw = div((Id- A)u) in fix(0,T), 

u = on 9fi x (0,T), 

w|i=o = u o i n ^- 



(5.4) 



Then keeping (3.15), (3.16), (3.17), and Proposition 1 in mind, we get for some constant 

SUp ||«(£)||wi ( m + ||Wi, V 2 M, VP||i 2 (nx(0,T)) 
0<<<T 

< C(||Id- A\\ Loc(nx{0yT)) \\u t ,V 2 u,VP\\ L2{ nx(o,T)) 

+ \\VAVu,A t u\\ L2{nx{0X)) + Hwollwjcn))- ( 5 - 5 ) 

The ^^(fi x (0,T)) regularity of the velocity, coming from (5.5), is not sufficient to 
control the Lagrangian coordinates, namely the terms containing A in the right-hand side 
of (5.5), because VW 2 ' 2 {Sl x (0,T)) is not embedded in Li(0, T; L 0O (Q)). Hence, to close 
the estimates, higher regularity is needed. Differentiating (5.4) once with respect to time 
is the easiest way to achieve it, because it does not affect the irregular density which is 
time independent in the Lagrangian setting. We get 



i]u tt - vAu t + V u P t = 

-v(A - A u )u t + u(A u ) t u - (V u ) t P in fix (0, T), 

&\v u u t = — div A t u in fix (0,T), 

M t = on dVt x (0,T). 



(5.6) 



At this stage the question of the regularity of Ut\t=o arises. This information can be 
found out only from the equations. At time t — the Eulerian and Lagrangian coordinates 
coincide (that is A = Id), so the regularity of u t \t=o is just that of r]~ 1 (VP — uAu)\ t=0 . 
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However the regularity of VP|i=o is unknown, so we rather have to use the fact that 
differentiating (5.4) 2 with respect to t implies that 

r]Ut\t=o + VP| t= o = z/Aw , divu t \ t =o = -div (A t \ t=0 u ). (5.7) 

Note that A t | t=0 need not be trivial so in order to bound u t \t=o in L 2 (fi), we first have to 
remove its potential part. For that, we use the Bogovskh operator B (see Lemma 5) setting 

= £>[— div (A t \ t=0 uo)] so that div0 = —div (A t \ t=0 u ) in Q, = at dQ. 

Let us notice that, because 

A(t,y)-DX(t,y) = Id and DX\ t=0 = Id with X(t,y) = Id+ / u{r,y)dr, 

Jo 

we have A|t=o — —Du , hence A t \ t=0 u = — (u ■ Vu ) = —div (u <g> u ). 

Now, W^i^i) is an algebra if n = 2,3. Hence A t \ t=0 uo is in Wl(VL) and the function 
defined above is in H / 2 1 (f2) 3 and satisfies: 



\w^(n) < C|l M o||^|(n)- ( 5 - 

Therefore System (5.7) recasts in 

r)(ut\t=o ~ 4>) + VP| t= o = vAu -r](p in fi 
div (u t \ t=0 — 0) = in f2 
(u t \ t =o -<P)\aa = on dtt. 
Now, testing the first equation by (u t \t=o — <fi) we get: 



T]\Ut\t=o — 4>\ dx < rj \uAu — r](f)\ dx. (5.9) 

n Jn 

Thus, due to (5.8), we discover that u t \t=o is in L 2 (f2) and that 

IKIt=o|U 2 (f2) < Cu,r,(\\u \\wi(n) + ll M o||^2 (Q) ). (5.10) 

At this point, we would like to apply an energy method to (5.6). However, as div u u t 
may be nonzero, one cannot eliminate the term coming from P t (which is out of control). 
So we modify (5.6) 2 by introducing a vector-field £ so that 

div u £ = -div ( A t u) in fi, _ ^ 

£ = on dfi. ' 

We need £ to satisfy suitable estimates (in terms of the right-hand side) in Loo(0, T; L 2 (f2))fl 
I/ 2 (0, T; W^i^i)) and £ 4 to be bounded in L 2 (fi x (0,T)). This may be done by means of 
a Bogovskh type operator construction as in [12]. Here we shall define £ (treating t as 
parameter) according to Lemma 6 in the Appendix. 



3 In fact, the function <p is in W|(0) but we shall not take advantage of this in what follows. 
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Let us start with the bound in L^O, T; L 2 (J1)): we have 

H^IUootO.T^n)) < \\A t u\\ Loo (0,T;L 2 (n))- 

Therefore, using the fact that 



A t = fe( A;+i )(- i ) fe+1 (/ D v ud A ) -a^, ( 5 - 12 ) 

we get (remember (5.1)) 

||^IUoo(0,T;L 2 (f7)) < ||«(8> VwllLoo^TjLa^)) 

< ll M IUo O (0,T;L oo (n))||VM|| Loo ( ,T;L 2 (n))- (5.13) 

In order to bound the right-hand side of (5.13), we apply the following classical parabolic 
estimate (which is related to our definition of Besov spaces in (1.3)): 



u 



lioc(o,T;B^- 2 / p (n)) ^ C{ y\\ u ^\\ B l;^(si) + IK V 2 w|| Lp(nx(0 ,T))). (5.14) 



Now, owing to Sobolev embedding, it is clear that the left-hand side of (5.14) controls the 
Lqo norm whenever 2 — Ijp > n/p, that is p > (n + 2)/2. The constant in (5.14) is time 
independent. Therefore, for any m G ((n + 2)/2,n*), we have 

IKK(£}x(0,T)) < C rn(\\u \\ B 2-2/ m(n) + \\u t , ^ 2 u\\ Lm (nx(0,T))) 

< C m (\\u \\ B 2-2/ n * (n) +T™~~S( UiP )(T)). (5.15) 

n* ,n* ^ ' 

Inequality (5.14) with p = 2 also yields 

\\Vu\\ Loo (p t T;L 2 (n) < IKIIw^cn) + IK, V u\\l 2 (Hx(o,t)) 

% ho\\w.^n) + T^^E {UtP) {T). (5.16) 

So, putting (5.13), (5.15) and (5.16) together, we get for some 5 > 0, 

UWL^T-Mm ^ IKIIw 2 2 (f2) + T -(u,P)( T )- ( 5 - 17 ) 

Next, in order to bound £ in L 2 (0,T; W^fi)), we use the fact that 

div u £ = -A t : Dm. 

In effect, owing to (3.13), one may write div (Au t ) = A : Du t , hence the above relation 
may be obtained by taking the time derivative of div (Au) = A : Du. So, using Lemma 6 
and remembering that n* > n, we get, for some 8 > 0, 

U\\L 2 (0,T;W^(n)) < C|||VM| 2 ||L 2 (f7x(0,T)) 

< C'||Vu||L oo (0,T;L 2 (n))||Vw||L 2 (0,T;I, oo (n)) 

< C{\\Vu \\ L2{n) +T 1 / 2 \\Vu t \\ L2{nx{0tT)) )\\Vu\\ L2 ^ T . w i t(n)) 

< CT s E (u , P) (T)(\\u \\ w i {a) +T s E (u , P) (T)). 
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Finally, let us bound £ t in L 2 (fi x (0, T)). For that, we apply the last part of Lemma 6 
which yields 

||6lU 2 (nx(o,T)) < C\\A t £i A tt u, A t ut\\L 2 ((ix(o,T))- (5.18) 

Using (5.13), (5.15) and (5.17), we get 

||A£|U 2 (nx(0,T)) < C'||VM|| L2 (o,T;L 00 (f7))|K||L 00 (0,T;L 2 (f7)) 

< C(ll«oll^(0) + T5 %, P) (T)) ||V U |U 2(0 ,T;^(n)), 

\\AtU\\L 2 (nx(0,T)) < ||Vu t ||i, 2 (nx(0,T))||w|U oo (0,T;L oo (f2)) 

< C||Vtit || L 2 (f2x(0,T))( || WO || ^2-2/"*^) +T E( Uj p)(T)), 

||Awt|U 2 (f2x(0,T)) < C\\Vu<S>U t \\L 2 (nx(0,T)) 

< C||Vw|| L2 (0,T;L oo (n))(||Wi - C|Uoo(0,T;Z, 2 (f2)) + ||ClUoo(0,T;L 2 (n)))- 

So one may conclude that 

||6||l 2 (Qx(0,T)) < (ll M o|| W 2 (Q) +T S (UiP) (T))||VM|| L2 ( 0iT . Loo ( n )) + Ulioll^a-a/n*^) 

+ T s E {U:P) (T))\\Vu t \\ L2{0tT . L2m + ||V«||L 2 (o,r;L 00 (ft))||t*t-^||L 00 (o,r;L 2 (n))- (5.19) 

Note that in (5.19), there is no factor T s for the leading order terms ||w t — ^||l 00 (o,t ; l 2 (q)) 
and ||Vw t ||z, 2 (nx(o,T))- 

Once the vector-field £ has been constructed, one may recast System (5.6) in 

V(ut ~ 0* - " A (ut - + V U P* 

= -z/(A - A u )u t + i/(A„) t u - (V„) t P - vtt + ^A£ in fix (0, T), 

div„(«t-O = in fix(0,T), (5.20) 

u t -f = on <9fi x (0,T), 

K~Olt=oe L 2 (fi) in ft. 

Note that, now, div u (u t — £) = and that (it t — £)l*=o is m L 2 (fi) with 

UK - Olt=olU 2 (ft) < C'(ll«o|lw^(n) + \\ u o\\w2 {n) ). (5.21) 

So taking the L 2 (ft)-inner product of (5.20) i with u t — £, there is no term generated by 
V u P t and we thus get 



T r 1 r 5 



\ I r]\ut-i?dx + u \V{u t -i)\ 2 dxdt<\ [ri\u t -£\ 2 dx +V/, 

^ </n *= T Jo Jn z Jn t=0 • 1 
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with 

h ■= Cv [ [ \(A-ld)\\Vu t \\V(u t -0\dxdt, 
Jo Jn 

h := Cv I I \Vu\ 2 \V(u t -£)\dxdt, 
Jo Jn 

I 3 ■= C I I \Vu\\VP\\u t -{,\dxdt, 
Jo Jn 

h ■= / r)\€f(u t -€)\dxdt, 
Jo Jn 

h ■= v \ I |V£- V(u t -£)\dxdt. 
Jo Jn 

In order to bound terms ii, J 2 , I 3 and I 5 , we use Holder and Young inequalities. We 

get for all e > 0, 

h < ev\\V(u t -OllL 2 (nx(o,r)) + C e ,u\\A - Id|| Loo ( nx ( 0iT ))||V-u t || L2(nx(0iT)) , (5.22) 

h < ^\\^(u t - ^\\l 2{nx{0tT)) + c £jU \\\Vu\ 2 \\l 2{nx{0tT)) , (5.23) 

h < £\\Ut - ^Wl^^T-M^)) +^H^ 7 ' U llL2(0,T;L oo (Q))ll^ 7 - P llL2(nx(0,T))' ( 5 -24) 

h <^||V(w t -Olli 2 (nx(o,T)) + C^I|VC||i 2 (n X (o,T))- ( 5 - 25 ) 

Inequality (5.24) deserves a remark : in order to "close the estimates" , we have to factor 
out the last term in the right-hand side by a quantity which is small enough when T goes 
to 0. Here this follows from the embedding W n l n , C L 2 (0,T; L^Vt)) which gives, because 
n* > 2, 

llVullLap.T^n)) < CT*-£\\u\\ w 2,i , ( nx(o,r))- ( 5 - 26 ) 



1 i 1 „„/\ X .— 2 



Finally, taking m G (1,2) so that 1 = — + -\, and S :— — — 1, we may write 

^4 < ||Wt-C|U n *(^x(0,T))||6|k m (f2x(0,T)), 

< e\\Ut — <H\L n *(ax(0,T)) + ^e^ llst|lL 2 (nx(0,T))- 

Combining interpolation and Sobolev embedding, we may write for all p G (2, oo), 

1 1 t*t - 4 \\L p (o,T;L q (n)) £ 1 1 tit - 4 \\ Loo (o,T;L 2 (n)) 11-^* _ ?J IK(fix(o,r))> 
with n/g = n/2 — 2/p. So taking p = q = n* := 2(n + 2)/n, we get 

IK - £lU„.(nx(0,T)) < C(IK - ClUoo(0,T;L 2 (fi)) + || VK - OIUa(nx(o,T)))- (5.27) 

Therefore, the above estimates for I\ to 1$ (with £ small enough) eventually imply that 

IK - £llL(o,r;L 2 (n)) + HI V K - Olli 2 (fix(o,T)) < 2 IIK - Olt=o||i 2 (n) 

+ C {j\^ U \\L 2 (0,T-,L ac (n)){\\^ u \\L oo (0,T-,L 2 (n)) + \\V P \\L 2 (nx(0,T))) 

1 1 . 1 1 O ii 1 1 o r 1 1 1 1 O ii ^_^ 1 1 O \ 

+ ll^ M l| Ll ( 0iT . Loo ( n ))||V-Ui|| L2 ( Qx ( 0iT )) +T ||?t||i 2 (fix(o,T)) + ll^slli 2 (nx(o,T)) )) 
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whence, using also the estimates for £, £ t and for Vw in L 2 (0,T; L 00 (fi)), we end up with 
IKIUooCo.t^W) + l|Vu t ||i, 2 (f}x(o,T)) < C«„(l + 2" S( UiP) (T)) . (5.28) 

Here 5 > 0. Let us also stress that C„ depends only on ||uo||w|(n)> IIpoIIl^), \\Pq 1 \\l x (q), 
Q and v. In particular, it is time-independent. 

Remark 5 At this stage, we find a limitation on the dimension of the domain: as we 
need to have Vw G L2{0,T\L oo (VL)), embedding requires that n* > n. This is fulfilled if 
n = 2 (because 2* = 4) or n = 3 (because 3* = 10/3,) but this is no longer satisfied 
in higher dimension. So we see that our method cannot be directly applied for n > 4 
unless we differentiate the system with respect to time, more times. Physical motivation 
for considering dimension n > 4 is unclear, though. 

Keeping (5.27) in mind, we see that in order to close the estimates, it suffices to bound 
the terms ||V 2 M||L n »(Q X (o,T)) and ||VP||z, n »(n x (o,T)) which appear in the right-hand side of 
(5.28). For that, we rewrite System (5.4) as a stationary Stokes system, treating r\u t as a 
source term, and the time variable as a parameter. So we consider 

-vAu + VP = -rju t - v(A - A u )u + (V - V U )P in fix(0,T), 

divw = div((Id- A)u) = -A : Du in fix(0,T), (5.29) 

u = on dVt x (0,T). 

Note that one may use Proposition 3 so as to handle the potential part of u. Therefore 
using standard results for the stationary Stokes equation (see [15]) enables us to get 

||zA7 2 w, VP||z, n «(n X (o,T)) < C(||ii t ||L f ,.(nx(o,r)) 

+ \\v(A - A>, (V - V U )P, vV{A : Dw)|| Lfi . (nx(0l T))). (5.30) 

The key to bounding the right-hand side is that, because n* > n, we have by embedding 
and Holder's inequality 



f Du(t',y)dt'\\ Lo 
Jo 



(Qx(0,T)) < CT l ^IIVuIIl^o.T;!,^)) 

< CT ~^\\ V u||x,„,(nx(o,T)) ■ 



In particular, this allows to write that 

|| (A - A u )w|| Ln „(nx(0,T)) < ||P(-4 T A)||L oo (0,T;L„»(Q)|P«|U n *(0,T;L oo (n) 

+ \\Id-A T A\\r ( n x( r ) . T)] \\D 2 u\ 



\L oo (Clx(0,T))\\^ f «||L„*(f2x(0,T)); 



< IH> 2 «||Li(o,r i i,„.(n)|pu||L n ,(o 1 r;L 00 (n)) 



rs^ 



\Du\\L,(nT:L~(n-\\\\D 2 u\ 



\L 1 (0,T;L oo (Cl))\\^' U\\L n »(Clx(0,T)), 

Similar estimates hold true for the other terms of the right-hand side of (5.30). So finally, 
putting together all the above inequalities leads to 

S (U ,P)(T) < C U0 (1 + T 5 S (UiP) (T)) 3 (5.31) 

for some 5 = S(n) G (0, 1) which may be computed explicitly. 
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Then we are able to close the estimate, namely to write that 

S (til P)(T) < 8C U0 (5.32) 

whenever T has been chosen so that 

8C U() T 5 < 1. (5.33) 

5.2 The proof of existence 

In this short subsection, we explain how the proof of existence may be achieved from 
the above a priori estimates. 

Taking for granted the proof of the existence of a solution in the smooth case is the 
shortest way. Under the assumptions of Theorem 2, one may for instance smooth out the 
initial density po by convolution by a positive mollifier. This provides us with a family 
of smooth approximate densities (po)e>o satisfying the same lower and upper bound as 
po- Then applying the local and existence and uniqueness statement of e.g. [6] (bounded 
domain case) or [8] (whole space case), one obtains a family of solutions (p £ ,v £ , VQ £ ) for 
System 0.1 with data (pq,vq). This family of solutions has the required regularity. In 
addition, the possible blow-up of (p £ , v £ , VQ £ ) at time T is controlled by the norm of Dv £ 
in Loo(0,T; L 2 (£l)) l~l Li(0,T; L 00 (Q)). Note that Proposition 2 ensures that (p £ ,v £ , VQ £ ) 
corresponds to a solution (i] £ ,u £ ,'VP £ ) of (3.10) with the same regularity. 

Now, the computations that have been performed in the previous section, combined 
with the aforementioned blow-up criterion ensure that the lifespan of (i] £ ,u £ , VP £ ) (or of 
(p £ ,v £ , VQ £ )) may be bounded by below as in (5.33), and that (5.32) is satisfied. The 
important point is that all those bounds depend on the density only through its infimum 
and supremum. So eventuality, (p £ ,v £ , VQ £ ) is uniformly bounded in 

Z,oo(n x (0,T)) x W$ n .(to x (0,T)) x L n *(Q x (0,T)), 

and in addition, d t v £ is bounded in Loo(0,T; L 2 (J1)) fl L 2 (0,T; W^ift)). 

By resorting to standard compactness argument, it is now easy to conclude that this 
family converges, up to extraction, to some (p, v, VQ) with the same regularity and satis- 
fying the same bounds. The regularity is so high that the it is clear that it satisfies (0.1). 
Uniqueness then follows from Theorem 1. 

Remark 6 An alternative approach to the issue of existence can be done by an iterative 
scheme performed in the same way as in our recent work [11], or as in [27, 31] for the 
homogeneous Navier-Stokes equations in the Lagrangian coordinates. 

6 Global existence 

This section is dedicated to the proof of global-in-time solutions. As pointed out in the 
introduction, in the case of smooth data, this is a classical issue that has been solved by 
different authors in the Eulerian framework : if there is no vacuum initially then global 
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existence may be achieved for general (smooth) data in the two-dimensional case and if 
the velocity is small in the three-dimensional case (see e.g. [8, 21]). 

As for us, in order to show the global existence, one may adopt the Eulerian approach, 
too. However the very low regularity of the density will enforce us to treat the inhomogene- 
ity of the fluid as a perturbation (hence to assume (1.7) or (1.9)) and to use the Lagrangian 
framework to prove the uniqueness. 

6.1 The two-dimensional case 

Here we prove Theorem 3. We concentrate on the proof of global a priori estimates. 
Indeed, existence can be established by an elementary approximation with smooth density 
exactly as in the previous section : for smooth enough densities, the existence of global 
solutions with velocity (locally) in W A ' 2 (Q x R + ) is ensured by [7] (bounded case) or by [8] 
(whole space case). In addition, let us emphasize that, in dimension two, this regularity 
guarantees that we are allowed to change coordinates between the Eulerian and Lagrangian 
ones. So the uniqueness follows from Theorem 1. 

For getting the global existence, the computations are simpler in the Eulerian frame- 
work. We aim at getting a control over vt in Loo(0, T; 1/2(0)) in terms of the data and of T 
only. Even though it is classical (see e.g. [3, 8, 21]) we here recall how to proceed. First, 
we test the momentum equation of System (0.1) by v t . We get: 

/^l 2 <fa + 5#/|v»|'«fa+/^-(V?»-v«)<b = o. 

Jn 2 at J n J Q 

Hence Holder and Young inequalities imply that 

IIVMHLw + ^II^IlL^) < IIVHlLwIlVvlll^n)- (6.1) 

On the other hand, using maximal regularity for the stationary Stokes equation 

— uAv + VQ = y/p( yfpv t + y/pv ■ Vi> ) in f2 

divf = in Q 

v = on dfl, 

gives (omitting the time-dependency) 

H|V 2 v|| L2( n) + ||VQ|| L2( n) < llv^lUooCn^llv^tllMn) + IIVHIm^H^Hm^))- ( 6 - 2 ) 
Now applying Ladyzhenskaya inequality ||Vi>|| 2 /^ < ||Vi> ||L 2 (^)||V 2 f ||L 2 (ft)> yields 

H|V 2 v||L 2 (n) + ||VQ||L a (n) < \\VPhoo(n)\\^/pvt\\L 2 (n) 

llplUooffi)llVHIi 4( n)l|V^lk(n) 
+ - [b.6) 

Making use of Ladyzhenskaya inequality in (6.1), also leads to 

d IIVV 2 



, \\\7 v \\ 1 1 1 1 

II l— I|2 , a ||V7 Il2 ^ 2 11 \\L2(£l) . ,-, \\P Loo(f^) .I /— ||4 Mv7 ||2 /r a\ 

\\Vpvt\\L2{(i) + ^llVvll^n) < v -r- + C \\Vpv\\ Lm \\Vv\\ L2{n) . (6.4) 

at \\p r„rni v 
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Finally, adding up (6.4) and (6.3), using that ||p(t)|U 00 (n) = HpoH^a) and performing a 
time integration yields 



2( } Jo V 2i/ HlPolk^) ||A)IUoo(n)/ 



^ II V7 112 , ^P kco(fi) /" II I- 114 || V7 112 J 



i/° 



o 



and Gronwall lemma implies that 



H|Vu(t)||i 2(n) + / (Hv^tllLcn) + llA)|IZi(ft)l|VQ||i 2( n) + ^IIPo|IZi ( n)llV^|lL(n))^ 

< HIVvoIImo) e C "~ 3||polU ~< n >-tf H^ll^cn)*-. (6.5) 

Note that the exponential term is controlled thanks to the basic energy equality (1.1) 
(combined with Ladyzhenskaya's inequality). 

Since W r 2 1 (IR 2 ) is not embedded into L^M 2 ) we still do not control the change of co- 
ordinates so that we cannot apply Theorem 1 to get uniqueness. So we are required to 
improve the regularity of the solution to (0.1). In fact, it turns out to be possible to obtain 
Wq'p smoothness for any 1 < p < oo and n < q < oo via bootstrap method. To avoid 
technicality, we focus on the case p = 2 and q = 4 which suffices both to perform the 
change of coordinates and to apply the uniqueness result stated in Theorem 1. We rewrite 
System (0.1) as 

mvt — uAv + VQ — (m — p)vt — pv ■ Vi> in Q x (0, T), 

divw = in fix(0,T), (6.6) 

v = in <9fix(0,T). 

where m = mf y£ Q p (y). Note that the method of characteristics ensures that the initial 
density controls lower and upper pointwise bounds of the density over Q x (0, T). 
Then using Theorem 5 we get: 



sup Vmv\\v(t)\\ B i mi) + \\mv t ,uV 2 v,VQ\\ L2{ o,T;L 4 (n)) 

0<t<T 

< C(\\(p - m)v t \\ L2 (o,T;L4n)) + \\pv •Vv\\ L2 ^ t -l a ^)) + V^lkollsi^a))- (6.7) 
Now, we have 

\\pv ■ Vw|| L2 (o,T;L 4 (n)) < C||u|| L4 (o,r;L 00 (n))||Vu||L 4 (n X (o ) r))- (6-8) 

The right-hand side of (6.7) is bounded by means of (6.5) as 

The first term of the right-hand side of (6.7) can be absorbed by the left-hand side provided 
c is sufficiently small in (1.7). This enables us to justify that the velocity v remains in 
W^(tt x (0,T)) for all T > 0. 

Finally, as 4 > 2 = dimf2, we are allowed to apply Theorem 1 in order to get the 
uniqueness of our constructed solutions. Theorem 3 is proved. 
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6.2 Global existence in the n-dimensional case 

In this part we address the global solvability issue in bounded n-dimensional domains 
with n > 3. We adopt the Lagrangian framework (however the Eulerian framework may 
be used as well, as regards the existence theory). The result is based on the technique for 
the homogeneous system performed in [26]. 

In contrast with the other sections, working in bounded domains is important: this is 
due to the following result which ensures the exponential decay of the energy norm. 

Lemma 4 Let u be a sufficiently smooth solution to (3.10). Then 

1 d f • 12 j f IV7 12 



-- / V \u\ 2 dy + u \V u u\ 2 dy = 0, (6.9) 

as long as the Lagrangian coordinates are defined. In addition if Q is bounded then 

IKOIlL^e-^ikllLw- (6-10) 

where Ai stands for the first eigenvalue of the Laplace operator, and rf = \\r]\\L 00 (a,)- 

Proof: The proof is similar to that of Lemma 1 : testing (3.10)2 by u we get (6.9). 
In order to get (6.10), it suffices to notice that, owing to incompressibility, we have 

r]\u\ 2 dy = / p\v\ 2 dx and / \V u u\ 2 dy= I \Vv\ 2 dx. 
Jn Jn Jn 

Hence using (1.1) and Poincare's inequality, we readily get (6.10). ■ 

Proof of Theorem 4: 

We focus on the proof of global a priori estimates for smooth solutions to (3.10). Indeed, 
from those estimates, it is easy to proceed as in Section 5 so as to prove the existence of 
a global solution under the assumptions of Theorem 4: this is only a matter of smoothing 
out the initial density so as to construct a sequence of smooth solutions (given by e.g. [7]) 
with uniform norms. 

So given a global solution (i],u, VP) to (3.10) with data p £ -^oo(^) satisfying (1.9) 
and u G B q ^, (£1) with divu = and u \qq = 0, we introduce the following quantities: 

M_i := ^ 1/p ^ 1/p 'lko|| B 2-2 /P(n) +m\\u \\ L2{nh 
M k := ^ 1/ ^ W ll«ll Loo(fcifc+1 . B ^/p (c)) + \\mu t ,uV 2 u,VP\\ Lp(k , k+1 . tLqm 

where m := inf p Q and k EN. Recall 1 < p < oo, n < q < oo. 
Let us notice that setting 

u(t,x) = vu(vt,x) and P(t,x) = v 2 P{yt,x) 

reduces our study to the case v — 1. Hence we shall assume from now on that v — 1. 

25 



Define a smooth function £ : R — > [0, 1] such that 

c« = { I * ; | \ (6.H) 

and set C fc (^) := C(* - k) for k > 0, and 4 := [A; - 1, A; + 1] for fc > 1. 
We recast System (3.10) with to = k — 1 as follows: 

m[( k u] t - A[( k u] + V[C fc P] = C fe (™ - ?7)«* 

-(A-A u )[( k u] + (V-V u )l( k P]+m(( k ) t u in fix(0,T), 

div [C fe w] = div {( k u] - div u [( k u] in fix (0, T), (6.12) 

( k u = on <9fi x (0,T), 

C fc M|j= fc _i = in fi. 

Let m* := suppo- We claim that there exist two positive constants K and a depending 
only on m*,n,Q,p,q, so that, under Condition (1.9), we have 

M k < KM_ ie ~ ak for all fceE (6.13) 

Let us observe that, by Sobolev embedding (here we use that q > n), we have 

/fc-i-i _ _ _ . 

\\Du\\ Loo{ n)d S < CJ2 \\D 2 u\\L p{h -,L qm <CY,M k . (6.14) 

£=0 £=0 

So, given that 

^ l-e~ Q 

fc>0 

if we assume that M_i is small enough -a condition which is equivalent to the smallness 
of d in (1.9)- then (5.1) is satisfied on [0, k + 1] if (6.13) is satisfied up to k. 

Proving (6.13) will be done by induction on k. The first step, k — 0, is clear. This 
is a direct consequence of Theorem 5 applied to (3.10) on the time interval [0,1], and of 
estimates for A. 



Let us now take for granted Inequality (6.13) up to k — 1. In order to prove it for k, 
we shall estimate (( k u, ( k P) on the interval I k . For that, one may resort once again to 
Theorem 5. First, we bound the right-hand side of (6.12) in L p (I k ; L q (Q)): we readily have 

\\C k ijn-vi)u t \\ Lv{Ik]Lqm) < (m* -m)\\u t \\ Lp{Ik . Lq{n)) , 

|| (A - A u )[( k u]\\ Lp{Ik . Lqm < \\A T A - Id\\ Loc{ nxi k )\\V 2 u\\ Lp{Ik;Lqm 

+ \\V(A T A)\\ Lx(Ik;Lqm \\Vu\\ Lp ( Ik;Loc (n)), 

II (V - V„)[C fe -P]||L p (/ fc; L 9 (n)) < \\A - Id\\ Loo (si x i k )\\VP\\ Lp (i k;Lq (n)), 

\\m(( k )t u \\L P (i k ;L q (n)) < rn*\\u\\ Lp{Ik , Lq(n)) . 

Let us notice that, by interpolation and because Q is bounded, we have for some 9 G (0, 1), 

|M|Lp(J fc ;L,(n)) < C \\ D2 A e L v (I k] L q (n))\\ U tL p ll k -L 2 (n)Y 
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Therefore, taking advantage of (6.10) and of the definition of M k _\ and of M k , we get for 
some [5 > (depending only on Q, p, q, m* and for all e G (0, 1), 

m*\\u\\ Lp{h . Lqm < e{M k ^ + M k ) + CM-i^- (6.15) 

Next, we have to bound the left-hand side of (6.12) 2 : we have 

\\div((Id-A)l( k u])\\ Lp{Ik . w i m < \\A-ld\\ Loo{nxIk) \\V 2 u\\ Lp{Ik . Lqm 

+l|V>l||L 00 (/ fc ;L,(n))||V«||L I ,(/ fc ;L 00 (n)), 

\\d t ((Id - A)( k u) \\ Lp (i k ;L q (n)) < \\A - Id\\ Loc{ nx(o,T))\\(( k u) t \\ Lp{Ik . iLqm 

+ \\At\\L p (I k ;L 00 (n))\\u\\L 00 (I k ;L q (n))- 

Let us look at the quantities depending on the matrix A. Recall that 

ft 

i-i 



A' 1 = Id + / Du(s) ds, 
Jo 

so taking advantage of (6.14) and of the hypothesis that follows, one may write that 

l|Id -^IUoo(nx(0,fe+l)) < 2||-DM||L 1 (0,fc+l;L oo (n)), 

and a similar inequality for Id — A T A. Likewise, we have 

k k 

\\DA\\ Loa(Ik , Lqm <Cj2\\D 2 u\\L q (i e ;L qm <CY,M e (6.16) 

e=o e=o 

and 

11^11^(4^(0)) < C\\Du\\ Lp(IkMn) < C(M k ^ + M k ). (6.17) 

So finally, putting together all the previous inequalities and applying Theorem 5 to (6.12), 
we end up with 

\\m(( k u t ),V 2 (C k u),V(C k P)\\ Lp{Ik , Lqm 



(* k \ 



At this point, it is clear that one has to take a = f3. Note also that if M_i and the 
oscillations of the density are small enough then, taking e small enough too, the above 
inequality implies, up to a change of C, 



- ak ), 



M k < C(cM fc _i + M_ ie - 
Now, using the induction hypothesis (6.13) for M k _i, we deduce that 

M k <KM_ ie - ak (— + - 
k ~ l \ K K 

Therefore, we see that if we take K = 1C and assume that c has been chosen so that 
cCe a < 1/2 then we get (6.13) for M k . 
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Note that our proof is not quite rigorous as we did use (6.13) at rank k in the above 
inequalities. To make the argument work, it is just a matter of replacing the interval Ij, 
with [k — 1,T]. By continuity of the norms with respect to time, it is clear that (6.13) at 
rang £ < k — 1 ensures that the desired inequality is satisfied on [k, T] for any T close 
enough to k. Then resorting to a standard bootstrap argument allows to conclude to the 
desired inequality for M&. This completes the proof of Theorem 4. 

7 Appendix 

Throughout this paper, we used repeatedly the following well-known result for the 
divergence equation (see e.g. [15] and the references therein): 

Lemma 5 Let Q be a bounded Lipschitz domain of M. n . There exists a linear operator B 
which is bounded from L q (Q) to W q (Q) for all q G (1, oo) and such that for any f G L q (Q) 
the vector-field u := B(f) satisfies 

div-u = / in Q and u\qq = on dQ. (7-1) 

This result may be proved by means of an explicit formula - the Bogovskii formula - that 
provides a solution to the above divergence equation in the case where Q is star-shaped. 
In our paper, we had to use a more elaborate version of the above lemma, namely the 
following statement that has been established in [12]: 

Proposition 3 Let Q be a C 2 bounded domain. There exists a linear operator B acting 
on couples (R, () with R : Q — > M. n and ( : dQ — > K. which is continuous from L q (Q; M. n ) x 

W q ' 1 ' q (dil, E) to L q (£l,W n ) for all q G (l,+oo) and such that u := B(R,() satisfies the 
generalized divergence equation: 

- f u-V(pdx = - f R-V(pdx+ f ((pda for all G C°° (H). (7.2) 

Jn Jn J an 

If in addition div R G L q (Q) and R ■ n = then u := B(R, 0) satisfies 

divu = divi? in il and u\qq = on dQ, 
and the following inequality holds true: 

\\u\\w}<si) < C||div-R|| L ,(n). (7.3) 

Furthermore, if we also have divi? G W q (Q) then u is in W q (Q) and we have 

||«||wj(ft) < C||divi?|| Wq i(n)- (7.4) 

We claim that this statement implies Lemma 2. Indeed, we set u := B(R,0). Then it 
is is clear that (2.4) holds true. Then differentiating u with respect to time yields 

u t = B(R t ,0). 

Hence applying the first part of the above statement yields (2.5). 
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In Section 5, owing to our use of Lagrangian coordinates, it was natural to extend 
Lemma 2 and Proposition 3 to the twisted divergence equation, namely 

div^ u = div R in Q and u = on dQ 

with div j^u := div (Am). In particular, we used the following statement: 

Lemma 6 Let Q be a C 2 bounded domain. Let A G L^VL; R raxra ) be such that det A = 1. 
There exists a positive constant c such that if 

||A-Id|| Loo( n)<c (7.5) 

then there exists a map B^ acting on couples (R, () with R : Q — >■ R n and ( : dQ — > R 
which is continuous from L g (f2;R n ) x W q ~ /9 (<9f2,R) to L q (£l,M. n ) for all q G (1, +oo) and 
such that u := Ba(R,C) satisfies the generalized twisted divergence equation: 

- I Au-V(pdx = - J R-V(pdx+ [ ((f)da for all (peC 00 ^). (7.6) 

Jn Jn Jan 

If in addition div R G L q (Q) and R- n — then u := Ba{R, 0) satisfies 

div a u = A : Du = div R in Q and u — on dQ, (7-7) 

and (7.3) with a constant independent of A. 

Finally, in the smooth case, if the data R and A depend on a parameter t in some 
interval o/R with R t in L q (Q) and (7.5) satisfied for almost all t, then u fulfills: 

IMk(ft) < C(\\A t u\\ Lq(n) + \\Rt\\ Lq (n)). (7.8) 

Proof: The proof follows from Proposition 3: we consider the linear operator T defined 
by 

T(0=£, where £ := B((Id- A)£ + R,{). (7.9) 

Note that this definition and the fact that det A = 1 imply that any fixed point of T 
satisfies (7.6) (or (7.7) in the smooth case with ( = 0). Next, under Condition (7.5) with c 
small enough, one may apply the Banach fixed point theorem to T so as to get a solution 
to our problem. The reader may refer to [12] for more details. 

Concerning the proof of inequality (7.8), it suffices to differentiate once the equality 

u = B((ld-A)u + R,0). 

We get 

u t = B(-A t u + (Id - A)u t + R t ,0). 

So it is a mere consequence of the first part of Proposition 3. ■ 
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